Abstract 



In order to explain the lack of carbon stars in the Galactic bulge, we have made 
a detailed study of thermal pulse - asymptotic giant branch stars by using a 
population synthesis code. The effects of the oxygen overabundance and the 
mass loss rate on the ratio of the number of carbon stars to that of oxygen stars 
in the Galactic bulge are discussed. We find that the oxygen overabundance 
which is about twice as large as that in the solar neighbourhood (close to the 
present observations) is insufficient to explain the rareness of carbon stars in 
the bulge. We suggest that the large mass loss rate may serve as a controlling 
factor in the ratio of the number of carbon stars to that of oxygen stars. 



0.1 Introduction 



Stars with initial masses in a range of O.QMq < M < SM© go through the 
asymptotic giant branch (AGB) phase at the end of their hfe. An AGB star 
has two nuclear burning shells: one is the burning helium surrounded a carbon- 
oxygen core, and the other is the burning hydrogen which is below a deep 
convective envelope. There exists an intershell region between helium burning 
shell and hydrogen burning shell. This region is rich in helium and carbon but 
has little neno and oxygen. [1] There are two special features in AGB evolution: 
thermal pulse (TP) and mass loss. [3] In each TP, the outer layers expand due 
to the high helium burning luminosity. The strong expansion has an influence 
on extinguishing the hydrogen shell. Then the deep convective envelope can 
penetrate into the intershell region and mix with the products of internal nucle- 
osynthesis. This mixing event is named the third dredge up (TDU). Following 
the dredge up, the hydrogen shell is re-ignited due to the star contraction. The 
star enters a phase of quiescent hydrogen burning known as the interpulse phase. 
In the current understanding of the stellar evolution theory, the carbon stars 
cannot be formed before the TP-AGB phase, the phase which occurs after the 
first TP to the end of AGB. After each TP, the surface abundances of TP-AGB 
stars may be modified by the TDU. 4 Following the dredge up, a certain quan- 
tity of intershell matter is brought up to the surface, and the surface abundance 
of ^^C becomes enriched. Under the repeated actions of TPs and TDUs, the 
stars can become carbon stars with ratio C jO > 1 in the envelopes if they are 
efficient enough. In recent decades, many studies have shown that the carbon 
stars in the Galactic bulge are very rare, compared with the numbers in the so- 
lar neighbours. [5ji3j The recent study [7] indicated that the stellar population of 
the Galactic bulge does not contain any normal carbon stars except one candi- 
date. However, the rareness of the bulge carbon stars is still a mystery. Feast[7] 
showed the lack of carbon stars in the bulge could not be explained simply by 
an age effect or a high metal abundance and suggested that the oxygen over- 
abundance maybe modified the ratio of the number of carbon stars to that of 
oxygen giants. 

In this paper, using the synthetic models and the different parameters of the 
initial oxygen abundance, we simulate TP-AGB stars in the Galactic bulge and 
discuss the effect of oxygen overabundance on the ratio of the number of carbon 
stars to that of oxygen TP-AGB stars (0-rich stars). In addition, we attempt 
to find the other controlling factors in the ratio of the number of carbon stars 
to that of 0-rich stars, for example, the mass loss. In Section 2, we present 
our assumptions and describe some details of the algorithm. In Section 3, we 
discuss the main results and the effects of the different parameters. In Section 
4, the main conclusions are given. 
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0.2 The Model 



0.2.1 Basic Parameters of the Monte Carlo Simulation 

For the population synthesis of the single stars, the main input parameters of 
the model are: (i) the initial mass function (IMF); (ii) the lower and the upper 
mass cut-offs Mi and for the initial mass function; (iii) the relative age of 
the single stellar population; (iv) the metallicity Z of the stars. [H [51 HD] 
Our primary mass distribution takes the IMF as[TT| 



where ^{m)dm is the probability with which a star has a mass (in solar units) 
between m and m + dm; tuq — O.IMq; ai = 0.29056; 02 = 0.15571. 

The lower and the upper mass cut-offs are M\ = Q.IMq and = I2OM0, 
respectively. If the initial mass is larger than 4.0Mq, the stars hardly become 
carbon stars because of the hot bottom burning which turns ^^C into ^^N. At 
the same time, the absence of the TDU prevents the carbon star from forming 
when the initial mass is less than 1.5Mq. Therefore, we lay a strong emphasis 
on the stars with the initial masses ranging from I.5M0 to A.OMq in this paper. 

The average value of [Fe/H] in the Galactic bulge is close to the mean value 
for the solar neighborhood. [12] In our work, we take the metallicity as Z=0.02 
for convenience. The maximum age of the single stellar population is 15 Gyr. 

The recent work has shown that the star formation rate in the Galactic bulge 
is 10 — IOOMq yr~^.[T3l[T4j We assume that the stars with total mass of about 
5OM0 are born in the bulge every year. 



The stellar evolution from the zero age main sequence to the first thermal pulse 
is dealt with using the rapid evolution code by Hurley et al.[Tn] After the first 
thermal pulse, we use a synthetic model for TP-AGB. The changes of the chem- 
ical abundances on the stellar surface during the giant branch (the first dredge 
up) and the early AGB (E-AGB) phase (the second dredge up) can be repre- 
sented by the simple fitting formulae in Refs.[Tni HZ]. In the following sections, 
the fitting formulae mainly come from Refs.[Tni HZ]. They are comprised of a 
Levenberg-Marquart gradient descent iterative ^^-minimization code. [IB] 

The initial abundances 

It has been known that oxygen in the bulge M giants is overabundant about 
twice larger than that in the solar neighborhood. [18] In our work, we take the 
initial abundances as those in Ref.[l9] for Z = 0.02 except oxygen abundance. 
To study the effect of the oxygen overabundance on the ratio of the number of 
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0.2.2 Synthesised TP-AGB Evolution 
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carbon stars to that of 0-rich stars, the oxygen overabundance is taken as the 
corresponding value of coefficient 9 in each case involved in our model. The 
details are shown in Table 1. The following data are the initial abundances 
shown by the mass fractions in our model: 
iH=0.68720, 4jje=o. 29280, 

i2C=2.92293x 10-3, i3c^4 jqsoOx 10-^ 

i4n=8.97864x lO--*, i5N=4.14000x 10-^, 

^^0=8.15085x10-3 X e, "0=3.87600x10-^ 
20Ne=2. 29390 xlO-^, 22]v^g^;^_45200x 10-^. 

The data for other elements are not shown here, and they do not affect our 
results. The enrichment of ^^O is equivalent to the enhancements of a elements 
which change the stellar opacity. Then, the enrichment of ^^O may have some 
influence on the stellar evolution. Salasnich et al[20] have shown that the en- 
hancements of a elements have a very small influence on the evolution of the low 
mass stars (M < ^Mq). The carbon stars originate from the stars with masses 
lower than AMq. In our paper, we neglect the influence of the enrichment of 
^^O on the stellar evolution. While the abundance of ^^O increases, we divide 

all metal element abundances by — 002 order to keep Z — 0.02, where 

1^0 = 8.15085 X 10-3. 



Core mass during the first thermal pulse 

Using the rapid evolution code by Hurley et al,[T5] we can obtain the stellar 
mass (MiTp) during the first thermal pulse. The core mass during the first 
thermal pulse, Mc^itp, is[l] 

Me,iTP = [-pi(Mo -P2) +Pz\f + (p4Afo +P5)(1 - /), (2) 

where / = (1 + ') ^ ^ and Afo is the initial mass in solar units. The 

coefficients pi, p2, ps, P4, ps, pe and are shown in Table 6 of Ref.[l]. 



Luminosity, radius and interpulse period 

We use the prescriptions by Izzard et al. [16j The luminosity is taken as the value 
calculated from Eq.(29) in Ref.[16]. We define radius RasL — AnaR^T^ff, where 
a is the Stefan-Boltzmann constant and Tcs is the effective temperature of the 
star. The radius is taken as the value calculated from Eq.(35) in Ref. ^16,. The 
interpulse period Tjp is 

logio(Tip/yr) = a28(Me/M0 - bzs) - - 10^^« + 0.15A2, (3) 

where the TDU efficiency A is defined in Section 2.2.4, and the coefficients are 

- -3.821, 
628 = 1.8926, 

C28 = -2.080 - 0.353Z + O.2OO(Me„v/Af0 + a~ 1.5), 
^28 = -0.626 - 7O.3O(MeaTp/M0 - 0{^M,/Mq), 

where Monv represents the envelope mass, a is the mixing length parameter and 
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equal to 1.75, C = log(Z/0.02), and AMc is the change in core mass defined as 



The minimum core mass for the TDU and the TDU efficiency 

The TDU can occur only in stars with masses above a certain core mass M™™. 
Groenewegen and de Jong took M™'" as a constant O.58M0.pT| Karakas et al[T] 
have found that Af""" depends on the stellar mass and the metallicity. They 
gave a fitting formula 

A'C"' = ai+ a2Mo + asM^ + a^M^ , (4) 

where coefficients ai, 02, 03 and are shown in Table 7 of Ref.[T] and Mq is 
the initial mass in solar units. According to the carbon star luminosity function 
in the Magellanic clouds, Marigo and Girardi[22J have considered that the pre- 
dicted Af™'" in Ref. ,1 is high. In this work, we take the value calculated from 
Eq.dlD as Af™". 

It should be recalled that M^™—Mc. itp if a stellar initial mass is Afinitiai > 
4Mq or Af ™ < A^ciTP • [D US 

The TDU efficiency is defined as A = ^amh'' ' ^^^'"^ AAfdrod is the mass 
dredged up to the stellar surface during the thermal pulse, and AAfjj is an 
increment of the core mass due to hydrogen burning during the preceding in- 
terpulse period. A is a very uncertain parameter. Karakas et al|lj showed a 
relation of A to be 

A(A^) = A,„ax[l - exp(-iV/iV,.)], (5) 

where A gradually increases towards an asymptotic A^ax with the increase of N 
(the progressive number of thermal pulsation) , and N,- is a constant determining 
how fast A reaches Amax- In our work, N,- is taken as the value calculated from 
Eq.(49) in Ref. [16], which can reproduce the results for Nr in Table 5 of Ref.[l]. 
Amax is given by (see Eq.(6) in Ref.[T]) 

_ bi + b2Mo + b^M^ 

where coefficients bi, 62, ^3 and 64 are shown in Table 8 of Ref. [T]. 
In this paper, we take A as the value calculated from Eq.([5]). 



Intershell abundances 

During the every thermal pulse, the dredged mass AAfdrcd incorporates into 
stellar envelope. According to the nucleosynthesis calculations by Boothroyd 
and Sackmann, [23j Marigo and Girardi|22j gave the fitting of the abundances 
of ^He, ^^C and ^^O in the intershell region which lies between helium burning 
shell and hydrogen burning shell: 
for AA^c <0.025A'/o, 

inter j^^4fjg) ^ q 95 _^ 400{AMc)^ - 20.0AAfc, 

inte^xl^^C) = 0.03 - 352(AAfc)^ + 17.6AA//c, 
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inter j^(-16q) ^ _32(AMc)2 + 1.6AMc, 

for AMc >0.025Mq, 

inter j^(4jjg) ^ Q 7Q ^_ o.65(AMc - 0.025), 

inter j^(12Qj ^ q 35 _ 0.65(AMc - 0.025), 
inter j^(160) = 0.02 - 0.065(AMc - 0.025). 

All the other isotopes are set to be zero in the intershell region and the abun- 
dances are renormalized such that their sum is 1.0. 



The third dredge up and the hot bottom burning 

If the hydrogen envelope of an AGB star is sufhciently massive, the hydrogen 
burning shell can extend to the bottom of the convective region. This process 
is named the hot bottom burning (HBB). For the HBB, we use a treatment 
similar to that in Ref. (21j. According to the model by Iben and Renzini.[24] 
Groenewegen and de Jong [21] gave the most suitable parameters for the fraction 
of newly dredged up matter exposed to high temperatures at the bottom of the 
envelope /hbb = 0.94, the fraction of the envelope matter mixed down with the 
bottom of the envelope /bur = 3 x 10^** and the exposure time of the matter 
in the region of the HBB tnBB — 0.0014rip. The temperature at the bottom of 
convective envelope Tbcc is given by Eq.(37) of Ref.|16j. 

For the drudged up masses, the quantities of materials added to the envelope 

are: 

A^He = '"'<=^A:(4He)AAfdrod, 
A12C = [(1 - /hbb)'"'°'^X(12C) 

+ tff/o™^™^Wd^]AMdred, 

A"C= [i^/J™-X™B(t)dt]AMd.ed, (7) 

A"N= [^/J™-X™B(i)di]AAfdred, 

A160 = [(1 - /hbb)'"'°'^X(160) 

+ fef/o™^™^Wd^]AAfdred, 

where X™^{t) is the chemical abundance of the material undergoing the HBB, 
and it is calculated in the way of Clayton's CNO bicycle. [25^ The details of 
Clayton's CNO bicycle can be seen in Ref.[25] (also see Refs.flHEJ). The CNO 
bicycle can be splited into the CN cycle and the ON cycle. The timescales in ON 
cycle are many thousands of times longer than those required to bring the CN 
cycle into equilibrium. Even in the most massive AGB stars undergoing vigorous 
HBB, the ON cycle never approaches equilibrium. Therefore, the effects of HBB 
mainly turn ^^C into ^''N and the abundance of ^®0 is not changed much. In the 
calculation of Clayton's CNO bicycle, the density on the base of the convective 
envelope is given by Eq.(42) of Ref. and the analytic expressions for the 
nuclear reaction rates in Clayton's CNO bicycle are cited from Ref . [26J . The 
initial conditions of X^^^{t) are X^^^{t = 0) = AT for and ^'^O, while 
= 0) = for "C and ^-^N. 
After every thermal pulse, the chemical abundances of stellar envelope X"''*s 
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each are 

^ X°'^M.„.(l-/,„.)+AX+ ^'7^^^— /;"^° Xit)dt (8) 

~ Monv+AAfdrod ' 

where the values of AX are given by expression ([7]) and the values of initial 
conditions of X{t) are the values of X{t = 0) that are equal to X°''*s. 

In this work, we care about the ratio of the number of carbon atoms to that 
of the oxygen atoms on the stellar surface and denote it as C/O. 



Mass loss 

The mass loss rate of the cool giant during AGB phase has a great influence on 
the chemical evolution of the stellar surface. In our work, two laws of the mass 
loss rate are considered: 

(i) A mass loss relation suggested by Vassiliadis and wood[27] based on the 
observations, and it is given as 

logAf= -11.4 + 0.0123(P ,s 
-100max(M/Mo - 2.5,0.0)), 



where P is the Mira pulsation period in days, given by 

(10) 



logP= -2.07+ 1.941og(i?/i?0) 



-O.9Olog(M/Af0). 
When P > 500 days, the steady super-wind phase is modeled by the law 

M(MQyr-i) = 2.06 X 10"^;^^, (11) 

where v^o is the terminal speed of the super- wind in km s^^; we use t;oo=15 km 
in this paper. 

(ii) A mass loss rate similar to Reimers' formula, which is given by Blocker pS] 
according to the simulations of shock-driven winds in the atmospheres of Mira- 
like stars l28j and expressed as 

M = 4.83 X 10-H'I-^-^L^-^Mji,^mers, (12) 

where MReimers is given by expression and r/ = 0.02. 30] 

In the other stellar evolutionary phases, the mass loss rates are given by 
Reimers' formula [5T| 

M = -4.0 X lO-i3^:^Moyr"\ (13) 

where L, R, M are the stellar luminosity, radius and mass in solar units, re- 
spectively, and the free parameter is ?7=0.5. 
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Table 1: Parameters of the models. 6 means the coefficient of oxygen overabun- 
dance. 



case 


coefficient 


Mass loss rate 


case 1 


1.0 


Eq.® 


case 2 


2.0 


Eq.® 


case 3 


3.0 


Eq.® 


case 4 


5.0 


Eq.® 


case 5 


2.0 


Eq.dni) 


case 6 


5.0 


Eq.iini) 



0.3 Results 

We construct a set of models with different input parameters. Table 1 lists all 
the cases considered in the present work. We calculate 1 x 10^ single stars for 
each case. 

0.3.1 The Effects of Parameters 

Positions in Hertzsprung-Russel (HR) -diagram In Fig.l, the different 
evolutionary phases during TP-AGB in HR-diagram are shown. For the same 
mass loss rate (see Figs.l(a)-l(d), or Figs. 1(e) and 1(f)), it is more difficult 
to form carbon stars with the increased of oxygen overabundance coefficient 6 
because the stars must drudge up more ^^C to the stellar surface in order to 
obtain C/0 > 1. When 9 — 5.0, the initial oxygen abundance is too high so 
that the dredged up ^^C is not enough for C/0 > 1, thus there is no carbon 
star (see Fig. 1(d)). On the other hand, for the same coefficient 6 (Figs. 1(b) 
and 1(e), or Figs. 1(d) and 1(f)), the mass loss rate taken as the value calculated 
from expression ® is more favorable for preventing the carbon stars for forming 
than that described by Eg. (fT2| . As we can see from Fig.l, the stars with the 
initial masses from 2.OM0 to A.OMq may become the carbon stars. 

Variety of C/O Figure 2 shows the varieties of C/O during the TP-AGB 
phase. For a star with an initial mass I.5M0, C/0 is almost a constant as 
shown in Fig. 2(a). Figures 2(b)-2(f) show, for the stars with initial masses 
ranging from 2.OM0 to 4.OM0, all of C/Os are changed but their varieties are 
different form each other. For a star with an initial mass 2.OM0 (see Fig. 2(b)), 
the core mass Mc is lower than Af™'" when it just evolves to the TP-AGB. It 
takes a long time to satisfy a certain core mass M™'" of the TDU. Therefore, it 
is hard to reach C/0 > 1. There is no carbon star with an initial mass 4.OM0 
(see Fig. 2(f)). The main reason is that the star with an initial mass higher than 
4.OM0 has a convective envelope which is thick enough, thus the temperature 
at the bottom of the convective envelope, Tbco, is high enough to undergo the 
HBB. The HBB turns ^^C in the drudge up materials into ^''N and prevents 
the carbon stars for forming. [iBj From Fig. 2(e), we can see obviously the largest 
variety of C/O in case 1 and the smallest in case 4 for the stars with the same 
initial masses. It means that case 4 is the most unfavorable case for forming the 
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Figure 1: The positions of TP-AGB stars in HR-diagram for the different eases. 
The sohd hnes mean 0-rich TP-AGB stars and the dotted hnes represent carbon 
stars. Stellar initial masses are shown by the figures in the side of the lines. 
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Figure 2: The varieties of C/O during the TP-AGB phase in the different cases 
and the initial masses. SteUar initial masses are shown in right top corners. The 
cases the different styles of lines correspond to are shown in the left top region 
of Fig.2(a). 



carbon star in our model. 

Lifetimes Figure 3 gives the stellar lifetimes of the different evolutionary 
phases. The phase of the carbon stars is a small fraction of the whole TP-AGB. 
From Fig. 3, we see that the peaks of lifetimes for the carbon stars each appear 
at 3.IM0, which are different from the peaks at 2.5Mq for the whole TP-AGB 
stars. The main reason is as follows: for a star with an initial mass 3M0, 
M™'" ~ Mc,iTP when it evolves into the TP-AGB phase. After a few thermal 
pulses the star begins to undergo the TDU. At the same time, according to 
expression JS]), the Amax is about 0.8 for a 3AIq star (also see Fig. 4 of Ref.[T]). 
Therefore, the TDU has a high efficiency A. The increases of its core mass and 
the mass loss of the envelope are slow so that the star takes a longer time to 
evolve to the end of TP-AGB. 

Figure 3 shows that the oxygen overabundance coefficient 6 has a significant 
influence on the lifetimes of the carbon stars. With the increase of coefficient 
0, the lifetimes of the carbon stars decrease. There is no carbon star when 9 
increases up to 5.0 (see Fig. 3(d)). In addition, Fig. 3 indicates that the mass 
loss rate also has an influence on the lifetimes of the carbon stars. For the same 
values of coefficient 9 (Figs. 3(b) and 3(e)), the lifetimes of the carbon stars with 
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Figure 3: The stellar lifetimes of the different evolutionary phases. The dot- 
dashed, dashed and solid lines represent the whole TP-AGB stars, 0-rich stars 
and Carbon stars, respectively. 



the mass loss rate as expressed by expression ^ are shorter than that given by 
expression (fT^ . 

Mass distribution Figure 4 shows the distribution of numbers (all the 
numbers are normalized to 1) of the modeled TP-AGB stars in the Galactic 
bulge as a function of initial stellar mass. Figure 4 shows that the initial masses 
of the progenitors of the carbon stars range mainly from 2.5Mq to 3.5Mq and 
the peaks of the initial masses of the carbon stars each appear at S.IMq, which 
is in good agreement with the case in Fig. 3. Comparing the cases in Fig. 4, we 
easily see the effects of the parameters 9 and the mass loss rate. The results are 
consistent with the descriptions in Fig. 3. 

0.3.2 Lifetimes and the Numbers of Carbon Stars 

In this subsection, we discuss the rough properties of the modeled population 
of carbon stars in the Galactic bulge and then proceed to make a more detailed 
comparison among the influences from the different assumptions. 

Table 2 gives the different models of TP-AGB population in the Galactic 
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Figure 4: The distributions of numbers (all the numbers are normalized to 1) of 
the modeled TP-AGB stars in the Galactic bulge as a function of initial stellar 
mass. The cases the different styles of lines correspond to are shown in the right 
top region of Fig. 4(a). 



Table 2: The different models of TP-AGB population. The first column gives case number 
according to Table 1. Columns 2, 4 and 6 give the numbers of TP-AGB stars, O-rich stars 
and carbon stars in the bulge, respectively. Columns 3, 5 and 7 show the average lifetimes 
of TP-AGB, O-rich and carbon stars, respectively. Column 8 is the ratios of the number of 
O-rich stars to that of carbon stars. 





TP-AGB Stars 


O-rich Stars 


Carbon Stars 


The ratios 


Cases 


Number 


Lifetime 


Number 


Lifetime 


Number 


Lifetime 


of 




10^ 


lO^years 


10^ 


lO^years 


10^ 


lO^years 


the number 


1 


2 


3 


4 


5 


6 


7 


8 


case 1 


113.23 


1.46 


107.03 


1.38 


6.20 


0.48 


17.3 


case 2 


113.23 


1.46 


110.94 


1.43 


2.28 


0.29 


48.7 


case 3 


113.23 


1.46 


112.59 


1.45 


0.63 


0.17 


178.7 


case 4 


113.23 


1.46 


113.23 


1.46 








00 


case 5 


138.52 


1.79 


135.07 


1.74 


3.45 


0.37 


39.2 


case 6 


138.52 


1.79 


138.51 


1.79 


0.02 


0.02 


6925.5 
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bulge. 

(i) The mass loss rate is taken as the value calculated from expression ([9]) 
(cases 1, 2, 3 and 4): the numbers of the carbon stars in the bulge may range 
from (case 4) to 6.02 x 10^ (case 1). Their average lifetimes range from 
(case 4) to 0.48 x 10^ years (easel). The ratios of the number of 0-rich stars to 
that of the carbon stars may range from 17.3 (case 1) to oo (case 4), where oo 
means that there is no carbon star in this case. The above results suggest that 
the oxygen overabundance has great effects on the numbers and the lifetimes 
of the carbon stars in the bulge. As shown in Table 2, the numbers and the 
lifetimes of the carbon stars decrease with coefficient 6 increasing. It is worth 
mentioning that Rich and Origlia[TB] indicated that in the bulge M giants oxygen 
overabundance is twice as big as that in the solar neighborhood. According to 
our model, — 2.0, which is insufficient for the lack of the carbon stars in the 
Galactic bulge. The carbon stars cannot be formed when 9 = 5.0 in our model, 
but the initial oxygen abundance is much higher than the observation. 

(ii) The mass loss rate is taken as the value calculated from expression l|12p 
(cases 5 and 6): Table 2 shows that the numbers and the average lifetimes of 
the carbon stars may range from 0.02 x 10^ (case 6) to 3.45 x 10^ (case 5) and 
from 0.02 x 10^ years (case 6) to 0.37 x 10^ years (case 5), respectively. It is 
obvious that the oxygen overabundance is favorable for preventing the carbon 
stars from forming, however, it is not enough to explain the lack of carbon stars 
in the bulge unless the oxygen overabundance coefficient 9 is equal to 5.0. This 
is consistent with the above result. 

(iii) For the same coefficient 9 (cases 2 and 5, or cases 4 and 6), we can 
see from Table 2 that the mass loss rate taken as the value calculated from 
expression ^ is helpful for explaining the lack of the carbon stars in the Galactic 
bulge. Recently, Willson[5^ gave a detailed description on d(logM)/d(logL) 
(see Fig.l of Ref.|32]). The average mass loss rate calculated from expression 
([9]) is higher than that calculated from expression ((T2|) during the TP-AGB 
in our simulations. The higher the mass loss rate is, the more quickly the 
envelope mass M^nv will decrease. According to expression the interpulse 
period Tjp increases with envelope mass decreasing. A long Tjp can reduce the 
TDU progressive number and TDU efficiency A. The large mass loss rate is 
unfavorable for forming the carbon stars. The environment in the Galactic 
bulge is significantly different from that in the solar neighborhood, for example, 
the higher irradiating luminosity and the stellar density. These factors may 
enhance the stellar mass loss rate. We suggest that the mass loss rate may be 
a factor of controlling the ratio of the number of 0-rich stars to the number of 
carbon stars in the Galactic bulge. 

0.4 Conclusion 

Employing the population synthesis, we make a detailed study of the TP-AGB 
stars in the Galactic bulge. We emphasize the relationship between the forma- 
tion of the carbon stars and the relevant parameters. The effects of the oxygen 
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overabundance cocfRcicnt and the mass loss rate on the ratio of the number 
of the carbon stars to the that of 0-rich stars are discussed. We find that 
the ratio of the number of the carbon stars to that of 0-rich stars is greatly 
affected by the oxygen overabundance coefficient 0. However, comparing with 
the present observations of the oxygen overabundance in the Galactic bulge, we 
believe that the oxygen overabundance which is about twice as large as that in 
the solar neighborhood is insufficient to explain the lack of the carbon stars in 
the bulge. In addition, we obtain that the mass loss rate also is a controlling 
factor in the ratio of the number of carbon stars to that of 0-rich stars. The 
larger mass loss rate is helpful for preventing the carbon stars from forming. 
There may be a more larger mass loss rate in the Galactic bulge than in the 
solar neighbourhood. Further work on this is desirable. 
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